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ABSTRACT 

The ages of stellar populations in 52 elliptical and SO galaxies in the Coma 
cluster are investigated using a new spectroscopic sample. More than 40% of the 
SOs are found to have experienced recent star formation in their central regions 
during the last ~ 5 Gyrs, while such activity is absent in the ellipticals. Galaxies 
in this sample have absolute magnitudes in the range —20.5 < Mb < —17.5, 
and the fraction of SO galaxies with recent star formation is higher at fainter 
luminosities. The observed luminosity range of SO galaxies with signs of recent 
star formation activity is consistent with them being the descendants of typical 
star-forming spirals at intermediate redshift whose star formation has been halted 
as a consequence of the dense environment. 

Subject headings: galaxies: clusters — galaxies: clusters: individual(Coma)— 
galaxies: elliptical and lenticular — galaxies: evolution 



1 Based on observations made with the William Herschel Telescope operated on the island of La Palma 
by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de 
Astrofisica de Canarias. 
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1. Introduction 

The formation mode and epoch of spheroids are among the most important unsolved 
issues in observational cosmology. These questions involve two possibly distinct aspects: 
when did the stars in spheroids form (i.e. their star formation history), and when were these 
stars assembled in a galactic structure that can be classified as a spheroid? 

The bulk of stars in spheroid-dominated (early-type) galaxies are believed to be old, at 
least in clusters (Bower, Lucey and Ellis 1992; Ellis et al. 1997; Kelson et al. 1997; Kodama 
et al. 1998; Stanford, Eisenhardt and Dickinson 1998; Kelson et al. 2001), but the frequency 
and amount of later star formation activity is still a matter of debate. To discriminate among 
different formation scenarios, it is crucial to determine whether field spheroidal galaxies had 
more extended star formation histories than their cluster counterparts (for opposite views 
see Bernardi et al. (1998) and Menanteau et al. (1999, 2001), also compare Renzini 1999 vs. 
Gonzales 1993 and Trager et al. 2000). 

In most studies, early-type galaxies (ellipticals and SOs) are treated as one class of 
objects, due to the general similarity in their colors and their "vicinity" along the Hubble 
sequence. This is in spite of the fact that also remarkably different formation processes 
have been proposed for ellipticals and lenticulars. In the case of ellipticals, the most cited 
alternatives are formation at high redshift through dissipationless collapse (the so-called 
monolithic scenario, e.g. Larson 1974), or a merger of disk systems (in the framework of 
hierarchical galaxy formation models, e.g. Baugh et al. 1996, Kauffmann 1996). Concerning 
the formation of SO galaxies, it is a longstanding issue whether they are "primordial" objects 
(i.e. they formed as SOs), or whether they were converted from spirals which lost their gas 
supply either by galaxy-galaxy collisions, ram-pressure stripping or gas evaporation in the 
hot intracluster medium, galactic winds or loss of their gas-rich envelopes (Spitzer and Baade 
1951; Gunn and Gott 1972; Faber and Gallagher 1976; Cowie and Songaila 1977; Burstein 
1979; Dressier 1980a; Larson, Tinsley and Caldwell 1980). 

A spiral origin for SO galaxies seemed to be favored by the discovery of large numbers 
of blue (presumably spiral) galaxies in clusters at z > 0.2 (Butcher and Oemler 1978, 1984), 
coupled with the abundance of SOs and ellipticals in clusters in the local universe (Hubble 
and Humason 1931; Oemler 1974). A renewed interest in studying the stellar population ages 
of ellipticals and SO galaxies separately arose in recent years, when the high spatial resolution 
imaging achieved with the Hubble Space Telescope (HST) uncovered a strong morphological 
evolution taking place in rich clusters during the last few Gyrs. There is an overabundance 
of spirals in cluster cores at z=0.5 and the SO galaxies are proportionally (a factor 2 to 3) 
less abundant than in nearby clusters, while the fraction of ellipticals at z=0.5 is the same 
or larger than at z=0 (Dressier et al. 1997). The progression of this morphological evolution 
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between z=0.5 and now has been recently clarified by Fasano et al. (2000), who explored 
for the first time the intermediate-redshift range z=0.f-0.25. These results strongly suggest 
that a large number of cluster spirals observed in distant clusters have evolved into the SOs 
that dominate rich clusters today (see also Kodama and Smail 2001). 

Observationally, the existence of systematic differences between the ages of stellar popu- 
lations in ellipticals and lenticulars is still controversial. Neither photometry nor spectroscopy 
of early-type galaxies in high-redshift clusters has been able to reveal any statistically sig- 
nificant difference (Ellis et al. 1997; Jones, Smail and Couch 2000), and the same is true 
for other works on lower redshift clusters (Jorgensen 1997, 1999; Ziegler et al. 2001; Trager, 
Faber and Dressier 2001). 

However, other studies have found evidence for young and intermediate age stellar pop- 
ulations in a significant fraction of lenticulars (Caldwell 1983; Gregg 1989; Bothun and Gregg 
1990; Fisher, Franx and Illingworth 1996). Recently, two studies have uncovered a difference 
in the recent star formation histories of cluster SO and elliptical galaxies. A high-quality 
line strength analysis has shown that most of the (faint) lenticular galaxies in the Fornax 
cluster have lower luminosity-weighted ages than the (brighter) ellipticals (Kuntschner and 
Davies 1998; Kuntschner 2000). In Abell 2218 at z=0.17, high-precision optical and near-IR 
photometry has found that ellipticals at all magnitudes and luminous SOs are old and coeval, 
while the faintest SOs appear to have younger luminosity-weighted ages (Smail et al. 2001). 

Both of these studies find that the differences between Es and SOs mostly stand out 
at faint magnitudes. The luminosity range explored by the different studies could be the 
key to understanding the reason for the apparently contrasting results found so far (Smail 
et al. 2001), but the relative role played by mass and morphology (SO vs E) in driving the 
evolutionary history of cluster early-type galaxies is still unclear (Kuntschner and Davies 
1998; Ziegler et al. 2001). 

The central issue is whether recent star formation in a significant fraction of the lentic- 
ulars is a widespread phenomenon in clusters. Additionally, it would be important to con- 
clusively understand if and, eventually, why this effect is evident only in a certain luminosity 
range. 

To address these questions we present here a study of the spectroscopically-determined 
ages and metallicities of ellipticals and lenticulars in what is considered the local prototype 
of rich clusters, the Coma cluster. Our sample comprises 52 galaxies (19 ellipticals and 33 
SOs) covering a broad range in luminosity (Mb ~ —20.5 to Mb ~ —17.5). We will show that 
the populations of SOs and Es differ for the presence/lack of recent star formation activity in 
a significant fraction of these galaxies. We also discuss the luminosity dependence of these 
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findings and its possible origin. 

2. Observations and analysis 

This work is based on the spectroscopic survey of galaxies in the Coma cluster presented 
in Mobasher et al. (2001, Paper II). Multi-fiber spectra with a resolution 6-9 A and a signal- 
to-noise ~ 15 — 19 were obtained for a random subset of galaxies from the sample of Colless 
and Dunn (1996) 2 in two areas of ~ 1 x 1.5 Mpc towards the center and the South- West 
region of the cluster. A full description of the sample selection, the observations and data 
reduction can be found in Paper II, where it is shown that this is essentially a magnitude 
limited sample with no significant bias. 

In this paper we present the stellar population properties of the elliptical (E) and SO 
galaxies. Morphological classifications from Dressier (1980b) are available for 77 galaxies in 
the Mobasher et al. sample, of which 19 are Es and 33 are SOs. All except two of these 
galaxies are in the magnitude range 13 < R < 16 (—20.5 < Mb < —17.5). Intermediate 
types such as E/SO and SO/E (2 galaxies) and SBOs (3 galaxies) have been excluded from 
the present analysis. Dressier classified as SOs those galaxies with a clearly recognizable 
non-spheroidal (disk or lens) component. When viewed face-on, they display an intensity 
discontinuity between the bulge and the disk. Galaxies with a smooth radial profile, with 
no intensity discontinuities, were classified as ellipticals. Disk galaxies with a clear spiral or 
outer ring pattern were classified as spirals. 

Line indices of the Lick/IDS system and emission line equivalent widths were mea- 
sured from the spectra and compared with spectrophotometric models to derive luminosity 
weighted ages and metallicities as described in Poggianti et al. (2001, Paper III). The line in- 
dices used in this paper, as well as those of the whole spectroscopic sample, will be published 
in a later paper of the series. Great care was taken to ensure an accurate calibration onto the 
Lick/IDS line index system, correcting for the spectral resolution, the galaxy velocity disper- 
sion and any residual offset. The comparison with standard Lick stars yields an uncertainty 
of ±0.1 A in our calibration of the H(3 index (see Paper III). As we will show below, this 
is significantly smaller than the random errors on our FL3 measurements, hence the latter 
dominate the age uncertainty. Emission-line spectra were excluded from the analysis of the 
Lick indices, and no emission correction was applied to the index measurements. 

To derive the stellar population properties, we have used the Padova version of Worthey's 



2 Complemented by an unpublished list of rcdshifts kindly provided by M. Colless. 
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models (see Paper III), which includes an accurate treatment of the horizontal branch of 
low metallicity stars. The main differences between the Padova and the standard Worthey 
models in the H/3-Mg 2 diagram that is employed below can be summarized as follows: a) the 
standard version does not cover the range [Fe/H]< —0.225 for ages < 8 Gyr, and b) the H/3 
strengths at ages 5-12 Gyr are generally lower in the Padova version, translating into age 
differences of typically 2 Gyr. At ages less than 4 Gyr, the difference is ~ 0.5 Gyr. We note 
that, while absolute ages are highly uncertain in this type of models, relative ages are much 
less affected by the model uncertainties. The current sample mostly consists of intermediate- 
luminosity galaxies, which do not display a high [Mg/Fe] ratio as more luminous ellipticals 
do (Worthey et al. 1992; Trager et al. 2000; Kuntschner et al. 2001, see also Paper III), 
therefore no attempt has been made to vary the abundance ratios of the models. 

We point out that: a) the spectra refer to the central 2.7 arcsec (1.3 kpc) of each 
galaxy, hence recent star formation episodes in the outer regions would not be detected; 
b) the quantities derived from the spectra are luminosity weighted ages and metallicities. 
Other quantities of interest - such as for example the mass fraction involved in the latest star 
formation episode - remain unknown; c) absolute ages/metallicities are much more uncertain 
than relative ages/metallicities, due to the uncertainties intrinsic to the models. We assume 
a distance modulus to Coma of 35.16 (< v >= 7000 kms -1 , H = GSkms" 1 Mpc -1 ). 

3. Results 

The main result is presented in Fig. 1, showing the H(3 versus Mg 2 diagram of ellipticals 
(empty symbols) and SO galaxies (filled symbols). The spectrum of one SO galaxy in our 
sample has very strong emission lines indicative of a current starburst and has been excluded 
from this and the following plots. 

There is a striking difference between the populations of ellipticals and SOs. All but one 
of the ellipticals are consistent within the errors with ages > 9 Gyr. In contrast, 13 out of 
32 SO galaxies have luminosity weighted ages smaller than 5 Gyr. 3 Hence, more than 40% of 



3 We note that only one of these - the galaxy with the luminosity- weighted age ~ 1.5 Gyr at H/3 = 3.15 — 
is classified as a "k+a" galaxy (="E+A") according to the MORPHs classification scheme by Dressier et al. 
(1999). Another 3 SOs with luminosity-weighted ages less than 5 Gyrs were classified as k/k+a's (marginal 
k+a cases). There is an excellent agreement between the MORPHs classification and the results based on 
the Lick indices for the strongest Balmer line cases. As expected, the MORPHs classification is able to 
identify those galaxies with luminosity-weighted ages < 2 Gyr regardless of metallicity but - due to the line 
strength thresholds adopted - it is insensitive to older episodes of star formation which, instead, can be still 
detected in the Lick system. 



- 6- 



the SOs appear to have experienced star formation at a recent epoch, while such widespread 
recent activity is absent in the elliptical population. These results are corroborated by the 
analysis of the H7 F - <Fe> diagram that provides a similar but independent method of 
deriving ages and metallicities. 4 

Inspection of Fig. 1 shows that the metallicity range covered by this sample is quite 
wide, typically from [Fe/H]=-0.7 to 0.4, i.e. abundances from 0.2 to 2.5 times solar. A 
metallicity spread is expected, given the broad magnitude range explored and the existence 
of a correlation between mean metallicity and galaxy luminosity. Ellipticals and SOs follow a 
broadly similar metallicity-luminosity relation, as shown in the top panel of Fig. 2. However, 
the relation for the ellipticals (Z = 1.850 — 0.137 x R, a = 0.180) is slightly flatter and less 
scattered than the relation for the SOs (Z = 5.721 - 0.400 x R, a = 0.284). This could 
explain the finding of Fisher, Franx and Illingworth (1996) that SOs might have a steeper 
Mg 2 -velocity dispersion relationship than Es. 

The bottom panel of Fig. 2 presents the luminosity-weighted ages as a function of 
magnitude. The population of SOs with young luminosity-weighted ages is clearly visible in 
the region of the plot below the dotted line. Only 4 out of the 13 SO galaxies with luminosity 
weighted ages < 5 Gyr are brighter than R = 14.5 (M B ~ —19). In general, the luminosity 
distributions of Es and SOs are different, with the SO distribution being more skewed towards 
fainter magnitudes. In fact, the E:S0 number ratio changes from 5:6 at 13 < R < 14 to 9:12 
at 14 < R < 15 and 5:13 at 15 < R < 16. 

It is therefore important to assess how luminosity effects come into play. This is better 
exemplified in Fig. 3, where the age distributions of Es and SOs are shown for the whole 
samples, and for galaxies brighter and fainter than R = 14.5 (Mb ~ —19) separately. 
The total age distributions of lenticulars and ellipticals differ significantly (with a 93.3% 
probability according to a Kolmogorov-Smirnov test). A higher proportion of "young" SOs 
in the faint subsample compared to the bright subsample is clearly visible in the middle panel 
of Fig. 3, and the difference between the age distributions of Es and SOs is more conspicuous 



4 During the course of this work, another spectroscopic survey of Coma galaxies has been presented 
(Castander et al. 2001). These authors investigate the star formation history of these galaxies using different 
methods, comparing with the spectral energy distributions of synthetic models. While a direct comparison of 
the "ages" derived would be meaningless - due to the different meaning of the word "age" between our and 
their analysis - it is interesting to note that in the great majority of cases the two studies agree in identifying 
galaxies with/without recent star formation: among the 21 galaxies in common (10 ellipticals and 11 SOs), 
both studies assign "old" ages (lack of recent star formation) to 15 galaxies, identify 1 emission-line galaxy 
in common and detect recent star formation (~ 3 Gyr or less) in 3 galaxies with absorption-line spectra. 
There is only 1 discrepant case of a galaxy in which we detect a recent activity not reported in their work, 
and another case where the results agree within the errors. 
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at faint magnitudes. 

This luminosity dependence of the evolutionary histories of SO galaxies is naturally ex- 
plained if a fraction of the SOs were previously star-forming spiral galaxies. We computed 
the luminosity evolution of a galaxy whose star formation was truncated at some time be- 
tween 2 and 5 Gyr ago. For galaxies with a star formation history (prior to truncation) 
typical of Sa, Sb, Sc and Sd galaxies (Barbaro and Poggianti 1997), we find that the B-band 
luminosity fades by 0.5 -1.5 magnitudes, depending on the truncation epoch and galaxy 
type. In the event of a starburst prior to the truncation, the luminosity evolution would be 
greater. Cluster spirals at intermediate redshift have M v ~ — 21 (Smail et al. 1997) in the 
cosmology adopted here, thus even in the most conservative scenario (only truncation, no 
starburst, B-V=0.4) 5 , the passive descendants of spirals are expected to be typically fainter 
than M B = —20/ — 19, corresponding to R = 13.5/14.5 in our Coma sample. 

Consequently, the great majority of galaxies that evolved from the spiral into the SO 
class should not be found at the bright end of the luminosity function. The luminous 
side of the SO luminosity function could be filled up by another population of SOs which 
did not evolve from spirals (or evolved from brighter spirals at higher redshifts). If the 
morphological transformation spiral =>- SO is accompanied by a conversion of a star-forming 
galaxy into a passive one, it is not surprising that those SOs with signs of recent star formation 
are preferentially lower luminosity galaxies. This could be the reason why many previous 
studies - concentrating primarily on the most luminous galaxies - have not found significant 
differences between the evolutionary histories of SO and elliptical galaxies. 

A question of great importance is whether it is mass, morphology, or environment 
that plays the dominant role in determining the evolutionary history of early-type galax- 
ies (Kuntschner and Davies 1998; Kuntschner 2000; Smail et al. 2001; Ziegler et al. 2001). 
The results presented here indicate that the distinction between an SO and an elliptical (i.e. 
the presence of a disk/lens according to Dressler's classification) is able to separate classes 
of galaxies with different star formation histories. These differences are noticeable between 
E and SO galaxies of similar luminosity. Hence, morphology appears to be the galactic 
property better correlated with the presence of recent episodes of star formation, at least 
in clusters like Coma. However, we have discussed how the differences between SOs and 
ellipticals can only be appreciated at Mb > —20, and become more and more evident going 
towards fainter magnitudes. In this respect, our findings for the Coma cluster are similar to 
those by Kuntschner and Davies (1998) for the Fornax cluster and by Smail et al. (2001) for 
Abell 2218. As far as galaxy metallicity is concerned, galaxy luminosity (i.e. mass) seems 



5 B — V ~ 0.4 for a very late-type spiral, ~ 0.9 for an early-type galaxy. 
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to be the primary parameter controlling it (see Paper III). 

While the effects of the cluster environment may be a viable mechanism to make spirals 
evolve into SOs, the existence of field lenticulars clearly proves that this cannot be the only 
mechanism of SO production. Two (or more) formation processes might be at work, of 
which one could be effective only in clusters and would be responsible for the formation of 
a fraction of the faint SOs. If this is the case and a large proportion of cluster SOs were 
spirals at higher redshift as suggested by the HST studies (Dressier et al. 1997), then the 
luminosity function of SO galaxies in distant clusters could be significantly different from 
the present-day luminosity function, likely with a lower number ratio of faint to bright SOs 
at high z. For the same reason, the field and the cluster SO luminosity functions at z=0 
would also be expected to differ. Studies of both the luminosity and the age distributions 
of SOs in different environments would be very valuable in placing constraints on acceptable 
evolutionary scenarios. 

In distant clusters, a large population of post-starburst/post-starforming galaxies in 
which star formation stopped some time during the previous 1.5 Gyr has long been known 
to be present (Dressier and Gunn 1983; Couch and Sharpies 1987; Dressier and Gunn 1992; 
Barger et al. 1998; Dressier et al. 1999, (MORPHs collaboration) and references therein). 
It has been suggested that these are recently infallen field galaxies whose star formation 
has been suppressed by the cluster environment. Furthermore, since most of these post- 
starforming galaxies have spiral morphologies at z=0.5, it has been proposed that at least 
some of these spirals are going to evolve into SOs at a later time, and that the timescale for 
morphological transformation (from spiral to SO) must be longer than the duration of the 
observational signature of the cessation of star formation (1.5 Gyr) (Poggianti et al. 1999). 
These post-starburst/post-starforming spirals would be an intermediate step between star- 
forming spirals and passive SOs. In this scenario, the fact that in Coma we find a number 
of SOs with "recent" star formation could appear to be at odds with the results found at 
z=0.5 where the post-starburst/post-starforming galaxies are mostly spirals. This apparent 
inconsistency is resolved if one considers that the Lick system employed here allows an 
exploration also of luminosity-weighted ages between 2-5 Gyr, so the times elapsed since 
the halting of the star formation that we can detect in this work are much longer than the 
post-starburst timescales which have been identified in distant clusters (1.5-2 Gyr at most). 
The SOs with recent star formation observed here could be a later evolutionary stage of the 
post-starburst /post-starforming spirals observed at z=0.5, but not as advanced as passive 
SOs observed in Coma and other clusters at z=0. 

Finally, it is interesting to examine the distribution in the sky of galaxies of different 
Hubble types and ages (Fig. 4). In our sample the phenomenon of recent star formation in 
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some lenticulars (filled symbols in the figure) is not preferentially observed in the South- West 
region around NGC 4839 6 . This seems to be in contrast with the findings of Caldwell et al. 
(1993); Caldwell and Rose (1997), who argued for a higher incidence of galaxies with recent 
star formation in the South- West area than towards the central region, and proposed that 
this phenomenon could be mostly related to the (possibly infalling) NGC 4839 group. Again, 
it is possible that the current study is sensitive to a wider luminosity-weighted age range 
than the studies by Caldwell and collaborators, and therefore identifies also star formation 
episodes that temporally preceded the last activity detected by these authors. The different 
epochs of activity could be associated with groups of galaxies in different regions of the 
cluster. 

Fig. 4 also shows that SO galaxies with luminosity-weighted ages < 5 Gyr have an 
asymmetric distribution with respect to the center of Coma, being preferentially located 
in a region east /north-east of NGC4874 (identified by the letter "A" in Fig. 4). Instead, 
older SO galaxies and ellipticals (triangles and circles with a central dot, respectively) are 
found to be spread out both north and south of NGC4874. This perception by eye of 
a difference in the spatial distribution of "young" and "old" SOs is not confirmed by a 
2D Kolmogorov-Smirnov test, which turns out to be inconclusive (the probability is only 
32.8%). Clearly, larger samples are needed for reaching statistically significant conclusions. 
The possibly higher level of dumpiness of the SOs with recent star formation could be due to 
the presence of substructure in Coma (Biviano et al. 1996; Colless and Dunn 1996; Zabludoff 
and Franx 1993, who found an apparently bi- or trimodal velocity distribution for Coma 
SOs), possibly being related to the NGC4889 merger event. At least some of the "young" SOs 
in the east /north-east region could be part of the same group, have accreted onto the cluster 
approximately at the same epoch, and consequently had their star formation halted by the 
effects of the cluster on similar timescales, as suggested by the fact they have comparable 
luminosity- weighted ages in Fig. 1. 

We are grateful to Matthew Colless, Ian Smail and Andrea Biviano for their comments 
and useful discussions regarding this work. 



6 See also Caldwell and Rose (1998) who found a similar result in their sample of low-luminosity early-type 
galaxies in Coma. 
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Fig. 1. — H/5 versus Mg 2 index strength for ellipticals (empty circles) and SO galaxies (filled 
circles). Overplotted are models from Paper III (see text). Random errors in the index 
measurements - taking into account the variance in each pixel and the statistical propagation 
of errors (see Paper III) - are shown. 



- 15 - 



-21 



-20 



-19 



-18 



-17 



M 



"5 n S3" 







E 

3 -1 



• • • 



o 
o 



-3 
30 



• o •••ooo mom m 



o 



20 



cn 1 

< 



5 

CD 



4 
3 



o 



o 



* • - ° • 



o 



o 



o 



* 8 . 

• o 



• m om 
m • 



12 



13 



14 



15 



16 



17 



Fig. 2. — Luminosity weighted metallicities (top) and ages (bottom) of ellipticals (empty 
circles) and SO galaxies (filled circles) as a function of the R-band magnitude. The absolute 
B magnitudes shown on top of the plot have been found for B-R=1.6. The observed R band 
magnitudes are aperture magnitudes over a radius 3 times the Kron radius and are taken 
from Komiyama et al. (2001, Paper I). Datapoints falling outside (on the right-hand side) 
of the model grid in Fig. 1 have been arbitrarily assigned a Z=3 in this plot. Similarly, 
an age = 30 Gyr was recorded when the datapoints lied below the model grid (see Paper 
III for a discussion of the mismatch between the observations and the model grid). The 
metallicity-luminosity relation of the top panel for all galaxies in this sample (Es+SOs) is 
Z = 4.555 - 0.318 x R, a = 0.238 (excluding the points lying at Z=3). 
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Fig. 3. — Luminosity- weighted age distributions of ellipticals and SO galaxies (left panel). 
In the middle and right panels the brightest (R < 14.5) and faintest (R > 14.5) subsets are 
plotted separately. 
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Fig. 4. — Positions in the sky of the early-type galaxies. Here we only show the area where 
there are galaxies in common between our and Dressler's sample. SO galaxies with luminosity- 
weighted ages < 5 Gyr are plotted as filled circles, while older SOs are represented by empty 
triangles. Empty circles are ellipticals. Small dots represent galaxies known to be members 
of the cluster from Colless and Dunn (1996). The positions of the two central dominant 
galaxies (NGC4874 and NGC4889) are indicated by the symbols "A" and "B", respectively. 
The cD galaxy in the South- West group (NGC4839) is identified with the symbol "C". 



